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ABSTRACT 

We present evidence for anomalous microwave emission (AME) in the Pleiades reflection nebula, us- 
ing data from the seven-year release of the Wilkinson Microwave Anisotropy Probe (WMAP) and from 
the COSMOSOMAS experiment. The flux integrated in a 1° radius around R.A.=56.24°, Dcc.=23.78° 
(J2000) is 2.15 ± 0.12 Jy at 22.8 GHz, where AME is dominant. COSMOSOMAS data show no sig- 
nificant emission, but allow to set upper limits of 0.94 and 1.58 Jy (99.7% C.L.) respectively at 10.9 
and 14.7 GHz, which are crucial to pin down the AME spectrum at these frequencies, and to discard 
any other emission mechanisms which could have an important contribution to the signal detected at 
22.8 GHz. We estimate the expected level of free-free emission from an extinction-corrected Ha tem- 
plate, while the thermal dust emission is characterized from infrared DIRBE data and extrapolated to 
microwave frequencies. When we deduct the contribution from these two components at 22.8 GHz the 
residual flux, associated with AME, is 2.12 ± 0.12 Jy (17.7er). The spectral energy distribution from 
10 to 60 GHz can be accurately fitted with a model of electric dipolc emission from small spinning 
dust grains distributed in two separated phases of molecular and atomic gas, respectively. 

The dust emissivity, calculated by correlating the 22.8 GHz data with 100 /jm data, is found to be 
4.36 ± 0.17 /iK/(MJy sr _1 ), a value considerably lower than in typical AME clouds, which present 
emissivities of ~ 20 /iK/(MJy sr" 1 ), although higher than the 0.2 /xK/(MJy sr _1 ) of the translucent 
cloud LDN 1780, where AME has recently been claimed. The physical properties of the Pleiades 
nebula, in particular its low extinction Ay ~ 0.4, indicate that this is indeed a much less opaque 
object than those were AME has usually been studied. This fact, together with the broad knowledge 
of the stellar content of this region, provides an excellent testbed for AME characterization in physical 
conditions different from those generally explored up to now. 

Subject headings: cosmic background radiation - diffuse radiation - ISM: individual objects (Pleiades) 
- radiation mechanisms: general - radio continuum: ISM 



1. INTRODUCTION 

Several microwave experi ments, like CO BE-DMK 
Kogut et al.|1996a|b|), OVRO ( | Leitch et al.|1997] Saska- 
toon (Td c Oliveira-C osta et al.|1997D, 19 GH z (de Ol iveira- 
Costa, et ai.||l998[) or Tenerife ( |de Oiiveira-Costa et al. 
1999), have revealed the presence of a statistical corre- 
lation between microwave maps and infrared maps trac- 
ing the thermal dust emission. These observations mo 
tivated a search for sim ilar signals in individual regions 
( Finkbeiner et aL] [2002), which led to unambiguous de- 
tections in several molecular c louds with different exper- 
iments, like COSMOSOMAS ([Watson et al.||2005l), AMI 
(|Ami Consortium: Scaifc et al. 2009a b), CB1 ( Casas- 
sus et al.||2006[ ICastellanos et al ||201ip, VSA Hl bbs et 



al| |2010| 



( jPlanck (Co llaboration et aT. 201ip 



or Plane 

These studies led to the idea that this dust-correlated 
signal, which was thenceforth referred to as "anomalous 
microwave emission" (AME), was indeed an additional 
diffuse foreground component, originated by a emission 
mechanism different from the well-known synchrotron, 
free-free and thermal dust emissions. 

AME, which is a significant contaminant of the cos- 
mic microwave background (CMB) in the range ~ 20 — 
60 GHz, was firs t thought to be free-free emission from 
T > 10 6 K gas (iLeitch et al.||l997|). This high temper- 



ature would be required to reconcile the observed mi- 
crowave intensity with the lack of intense Ha emission, 
as it would be expected in the case of bremsst rahlung 
emission with lower gas temperatures. However, |Draine| 
fc Lazarian| ( |1998a[ ) ruled out such high temperatures 



on energetic grounds, and proposed electric dipole emis- 
sion from small rotating dust grains in the interstellar 
medium (th e so-called spinn ing dust), a mechanism first 
proposed by 
for AME. T 



Erickson (1957), as a p lausible explanation 
re models provided by |Draine fc Lazarian] 
( 1998b ) for different media show peaked spectra with 



duce fairly well the observations ([Finkbeiner 


2004]_de 


Oliveira-Costa et al. 2004 Watson et al. |2005 


Iglesias- 


Groth 


2005, 2006; Casassus et al. 2006, 2008 I 


Dickinson 


et ai.|2009; Tibbs et al.|2010|). Recently, the previous an- 



introducing more detailed consideration s about the grain 
shapes and their rotational prop erties dAli-Hai'moud et 



al.||2009 Hoang et al. 



2010 



|Draine fc Lazarian 



explanation for AME 



2011 Silsbee et al. 



Q1999 ) presented an a 



2011) 



Iternative 

based on magnetic dipole radia- 
tion from hot ferromagnetic grains. Models of magnetic 
dipole emission based on single-domain grains, which 
predict polarization fractions significantly larger than 
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those of the electric dipole emission ( Lazarian fe Draine 
2000 ) , have been ruled out by different observations QBat 



tisteili et al."l20 06; Casa ssus et al.|20061|Kbgut et ai.|20U7 
Mason et al.||20091 [Lo pez Caraballo et al.||2011[ ). 

|Bennett et al.| (|2b03) proposed an alternative mecha- 
nism based on a fiat-spectrum synchrotron radiation as- 
sociated with star formation activity to explain Wilkin- 
son Microwave Anisotropy Probe (WMAP) first-year ob- 
servations. However, this hypothesis was questioned by 
other studies. De Oliveira-Costa at al. (2004) argued 
that the spinning dust hypothesis is clearly favored as it 



nal below 20 GHz 


Fernandez-Cerezo et al. 


2006 


Hilde- 


22.5 


brandt et al. 


2007 


, whereas Davies et ai. 


2006 


tound 



significant correlation m regions away t rom s tar-torming 
regions. More recently Kogut et al. (2011) concluded 
that spinning dust models can tit the ARCADE data at 
3, 8 and 10 GHz better than the flat-spectrum model. 

In this article we present a study of the AME aris- 
ing in the Pleiades reflection nebula, using data at ~ 
11-17 GHz from the COSMOSOMAS experiment and at 
^20 — 60 GHz from WMAP. This nearby region, having 
various bright stars capable of exciting the circumstel- 
lar medium, is an obvious AME candidate. Its relatively 
high Galactic latitude (b = —24°) makes it also attractive 
as we can largely avoid contamination from the Galactic 
plane. Furthermore, it has been widely studied, there 
being extensive information in the literature about its 
physical properties that allow us to better understand 
and to constrain the AME intensities. As we will see 
next, those studies normally refer to different dust sub- 
structures located within this region. However, the low 
angular resolution (~ 1°) of the microwave data we han- 
dle preclude us from studying those individual features, 
and therefore we will be focused in the whole complex, to 
which we will refer to as the "Pleiades reflection nebula" . 

2. PLEIADES ENVIROMENT AND PHYSICAL 
CONDITIONS 

The discovery of a reflection n ebulosity towards the 
Pleiades star cluster (Tempel 1861 ) revealed the presence 
of interstellar dust in this region. This dust is reflecting 
light from hot blue stars in this cluster, especially the 
bright stars 17 Tau, 20 Tau, 23 Tau and 25 T au, with B- 
band m agnitudes ranging from 2.81 to 4.11 (Ritchey et 
al.|2006[) and r espective spec tral types B6 III, al 111, B6 
IV and B7 III ( |White et aL]|2001| ). The Pleiades cluster 
lies within the Tau rus complex at a distance of 125 pc 
van Leeuwen||1999 ), and extends ~ 30' (~1 pc) around 

23°47' (J2000). |Caste- 



coordinates R.A.= 3 h 46 m , Dec 



laz et al.| ( 1987) discovered infrared dust emission around 
the positions of the Pleiades stars, which they attribute 
to nonequilibrium very small grains or to molecules with 
10-100 carbon atoms. The brightness distribution of this 
emission, which extends ~ 1°, and its connection with 
the exciting stars, can be appreciated in Figure Tl where 
we show a DSS blue-band image with 100 /im IRAS in- 
tensity contours overplotted. The strong emission ~ 15' 
south of 23 Tau is originated in the Mer o pe m olecular 



cloud, 
maps. 



identified by 
Another two 



Federman fc Willson| ([T984j) in CO 
aright knots ot infrared emission are 



seen in the IRAS maps concentrated around the stars 17 
Tau and 20 Tau. 
Except for the area of the Merope molecular cloud, 




Fig. 1. — Digitized sky survey blue-band photometrjQ in the 
region of the Pleiades star cluster with IRAS 100 fira data overlaid 
with contours. These contours correspond to intensity levels of 20, 
30, 50, fOO, 200 and 300 MJy sr _1 . Arrows point to the positions 
of the main exciting stars in the region, and the circle indicates the 
aperture we will use for flux integration. 

a http : //archive . stsci . edu/cgi-bin/dss_f orm This image was 
taken with the Oschin Schmidt 48-inch Telescope on Palomar 
Mountain. 



with .Eb-v = 0-35 mag measure d in the dire ction of 
the background star HD23512 by |White et al. 
the optical depth in this region is low (i/B_v ^ 
0.07 mag), and there s eem to be lit tle other dust along 
the cluster line of sight ( Cernis 1987 ) . The physical prop- 



(2001) 

umr- 



erties of the gas in different positions of th e whole com- 
plex have been studied by different authors. |Federman fc 
WiHson|(fi984[) obtained gas densities n ~ 300-500 cm~° 



and temperatures T g ~ 20 K from molecular line ob- 
servations in the Merope molecular cloud. They found 
column densities intermediate between diffuse and dark 
clouds, and concluded that Merop e is indeed m uch less 
opaque than typical Taurus clouds. White ( 1984 ) studied 
a larger area in the Pleiades enclosing the mam exciting 



stars and estimated d ensities 
~ 19841 



400 cm , whereas 
100 cm" 3 in a similar 



Gordon fc Arny| fll984|) fou nd n 
area. Ritchey et al.| ( |2006[ ) obtained even lower values of 
n ~ 46, 16 and 40 cm" J in the lines of sight of 25 Tau, 
27 Tau and 28 Tau, respectively. 
Different st udies, based on CO maps (Federman & 



Willson 1984) and on hig h-resolution spectra towards 
Pleiades stars ( White 2003 1 , have shown the gas and dust 
that make up the Pleiades nebula to have a significant 
radial velocity offset relative to the star cluster. A corol- 
lary of this is that the Pleiades stars were not formed out 
of this surrounding material but, more likely, the spatial 
association between the stars and the interstellar gas is 
the result of a chance encounter be tween the clus ter and 
one or more approaching clouds. White (1984) identi- 
fied two separated phases: one layer of atomic gas with 
(local standard of rest) radial velocity of 7 km s" 1 , ly- 
ing > 1 pc in front of the cluster and with thickness of 
< 0.3 pc, and one molecular phase (though strongly ion- 
ized) with radial velocity of 10 km s _1 within the cluster. 
The presence of these two phases was confirmed by latter 
studies, which also showed that the atomic phase proba- 
bly belongs to the Taurus clouds, but are inconclusive on 



Anomalous emission in the Pleiades 



3 



whether or not the molecular phase is embedded within 
the star cluster fWhite|2003] |Gibson fe Nordsieck||2003| ) . 



3. DATA ANALYSES 

3.1. Expected free-free emission from Ha 

The radio free- free emission must be accurately charac- 
terized in order to quantify any possible dust-correlated 
excess at microwave freq uencies. To this aim w e follow 



the method described in Dickinson et al. ( |2003 ), which 
allows to estimate the level ot this emission fr om Ha. We 
resort to the full-sky Ha map compiled by Fin kbeiner| 
(2003), a patch of which is shown in Figure [2] in the re- 
gion of the Pleiades reflection nebula. The solid circle 
depicts the 1° aperture we will use in this work for flux 
integration, while the area limited by the two dashed 
circles at radii 1.7° and 2.0° define the annulus where 
we will average the background signal that will be sub- 
tracted. However, in this case it is evident that this 
annulus would not give a reliable estimation of the back- 
ground, as it contains clear contaminant signal rather dif- 
ferent from the true background level around the Pleiades 
clouds. In fact, the averaged signal is larger in the back- 
ground annulus than in the aperture, giving a negative 
background-subtracted Ha intensity of —0.41 ± 0.09 R. 
When positioning the background annulus at different 
radii between 1.0° and 3.0° we consistently get negative 
values between « —0.2 R and sa —0.6 R. Therefore, from 
an aperture integration here we can only set upper limits 
on the Ha emission. Instead, we choose to proceed in a 
different way. 

Inside the aperture the only evident signal above the 
background level is a bright Ha knot located ~ 10' south 
of 23 Tau and coincident with the position of the Merope 
molecular cloud. On these grounds we consider that all 
the Ha emission associated with the Pleiades nebula is 
produced in this position. In order to calculate its associ- 
ated intensity we fit this feature with a two-dimensional 
elliptical Gaussian function plus a constant background 
level. We obtain a maximum amplitude (source plus 
background) of 5.90 R, a background level of 2.58 R, and 
FWHM of 14.3' and 20.7' in two orthogonal directions. 
Using these parameters the mean background-subtracted 
Ha intensity in the aperture is In a = 0.092 R. Note that 
this value is higher than the 3-sigma upper limit obtained 
through aperture photometry. Therefore, it is more con- 
servative, as it will be extrapolated to microwave frequen- 
cies in order to characterize any possible excess emission 
in that range. 
This value needs correction from Galactic extinction, 
Using 



rcorr . 



O'Donncll 



( 1994 ) poly- 



I Ha X 10^ H «)/ 2 -5. 

nomials that allow calculating the extinction at any given 
wavelength in the range ~ 0.3 — 0. 9 /mil, and taking 
i?v = -Av/^b-v — 3.6 from Cernis (1987), we obtain 
A(Ha)= 2.94 Eb-v- Ideally we should apply a pixel 
by pixel correction, for which we would need a redden- 
ing map. We have inspected the Schlegel et al. (1998) 
full-sky Eb-v map in this region, which gives a mean 
reddening inside the aperture of 0.30 mag, with a max- 
imum of 2.06 mag in a position ~ 15' south of 23 Tau, 
coi nciden t with the Merope molecular cloud. |Ritchey et] 
al. ( 2006 1 give reddening values towards twenty stars be- 
longing to the Pleiades (taken from White et al. 2001), 
with an average value of 0.074 mag, and a maximum of 



0.35 mag towards the star HD23512. These figures are 
consistent with other measur ements in the direction of 
Pleiades or backgroun d stars (White 1984 Cer nis||1987 



Herbig fc Simon | 2001[ ) . The reddenings from the |Schlegel| 
et al.| ( |1998j) map correspond to the full integral along 
the line ot sight, and therefore their higher values could 
be due to the extinction produced in background dust. 
However, given the relatively high Galactic latitude of 
the Pleiades, one would expect those differences to be 
smaller. Therefore, instead of performing a pixel by pixel 
correction, we conservatively assume an average value 
of -Eb-v = 0.1, which is slightly larger tha n the mean 
over the twenty stars of Ritchey et al. (2006). We finally 
obtain 1^" = 0.12 R, and assuming an electron tem- 
perature of 8000 K, we obtain an emission measure of 
EM = 0.267 cm -6 pc. From this, we derive the expected 
spectrum of free-free flux integrated over the aperture. 

3.2. Radio data 

There are different radio surveys from which we can 
also constrain the free-free emission. We use the IHaslaml 
et al.| ( |l982[ ) map at 0.408 G HzR the|Berkhuiisen|(|l972t 
map at 0.820 GHz and the |Reich fc Reich| fl986 \ map 
at 1.4 GHz[f] These maps smoothed to a common resolu- 
tion of 1° (except the 0.820 GHz map, whose original res- 
olution is 1.2°) in the position of the Pleiades are shown 
in Figure[2j No obvious excess above the background and 
noise levels is seen, so we have to establish upper limits 
for the flux densities at these frequencies. To this aim, 
we assume the noise per pixel to be equal to the standard 
deviation of the signal calculated in the background re- 
gion (external annulus) of the maps smoothed to 1° . The 
flux noises in the aperture derived from this assumption 
are respectively 0.38, 0.30 and 0.17 Jy (68% C.L.; in Ta- 
ble m we show upper limits at the 99.7% C.L.), while the 
free-free fluxes estimated from Ha (see previous section) 
at those frequencies are at a level of ~ 0.030 — 0.027 Jy. 

3.3. Microwave data 

3.3.1. COSMOSOMAS 

The COSMOSOMAS ( Cosmological Structures on 
Medium Angular Scales) experiment consisted of two 
circular-scanning instruments, so-called COSMOll and 
COSM015, operative between 1998 and 2008 from the 
Teide Observatory (altitude 2400 m, Tenerife). They 
produced 0.8-1.1° resolution daily maps with full cov- 
erage in right ascension and ^20° coverage in declina- 
tion in four frequency bands centered at 10.9 GHz for 
COSMOll, and 12.7, 14.7 and 16.3 GHz for COSM015. 
A description of the instrume ntal setup of COSMOll 
and COSMQ15 can be fou nd in |Hildebra ndt et al.|fl2007 1 
andlGallegos et al. (2001), respectively. 



An important step ot the data processing is the re- 
moval of the first seven harmonics of each scan, in order 
to suppress the 1// noise from the receiver and from 



1 We use the map supplied by 

2 Downloaded from 



Platania et al. 



(2003) 



http : //www.mpif r-bonn.mpg . de/ survey .html 

° These maps (and, in general, all used in this work) are pro- 
jected into the HEALPix ( [Gorski et al,|2005| pixelization scheme, 
using a pixel size of ~ 7 (N s ide = 012.J. In the case of the 
0.820 GHz data we carried out the re-pixelization ourselves, as 
they are given in a different format. 
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Fig. 2. — Ha intensity and radio temperature maps in the position of the Pleiades reflection nebula. Crosses mark the positions of the 
same stars shown in Figure [l] A small knot of Ha emission is seen ~ 10' south of 23 Tau, but in the radio maps its associated free-free 
emission is overshadowed by the noise and background fluctuations. 



the atmosphere. Each scan, which is a stack of 30 s of 
data (corresponding to 30 spins of the primary mirror) , is 
Fourier-fitted and terms up to seventh order (i.e. a con- 
stant term plus seven sin and cos terms) are removed. 
This results in a distortion of the broad angular scales 
and in a loss of flux. Further details on the data re- 
duction, calibration and map- making of COSMOll and 
COSM015 data can be found in 
and in 



Fernandez-Cerezo et al. 



Hildc brandt et al.| ( 2007 1 
1 2006|), respectively. In 



those references the final maps can be seen, as well as 
the discussion about the dust-correlated signal that was 
found at |6| > 20°. 

The maps we use here are a stack of ~ 150 days 
of observations in each of the two linear polarization 
modes of COSMOll (November 2003 to June 2005) 
and of - 110 days of COSM015 data (October 1999 
to January 2000). They are convolved to a common 
resolution of 1.12°, and have final sensitivities (aver- 
aged over the «10,000 square-degree observed sky re- 
gion) of 35 fjK beam" 1 for COSMOll, and 53, 56 and 
118 /iK beam -1 for the three frequencies of COSM015. 

COSMOll and COSM015 observed in different con- 
figurations, covering different declination ranges. The 
minimum declination was rs 23.5°. In the case of 
COSM015, there are also ~5 days in a configuration 
down to Dec. =16. 7°, but we discard these data because 
of being too noisy. Most of the area encompassed by the 
Pleiades reflection nebula is covered by the observations, 
but unfortunately its center lies just 17' from the border 
of the map. Also, the total number of days of obser- 
vations in the configurations reaching these low declina- 
tions are - 95 for COSMOll and - 12 for COSM015. 
However, despite the fewer observing days, the maps sen- 
sitivities in the region of the Pleiades are similar to the 
full-map average sensitivities quoted above thanks to the 
circular scanning strategy, which results in a better spa- 
tial coverage towards the edge of the map. Possible edge 
effects arisen in this region are avoided by ignoring the 
blank pixels in the convolution process. 

COSMOSOMAS maps in the region of the Pleiades are 
shown in Figure [3] No clear detection is seen in any of 
the frequency bands. Marginal evidence for a weak sig- 
nal in the position of the Pleiades could be claimed in 
the COSMOll map. We obtain an indirect estimate of 
0.67±0.09 Jy for the 11 GHz flux from the DIRBE fluxes 
calculated in section |3.4| by using the dust-correlated 



emissivities between the COSMOll map and the DIRBE 
100, 140 and 240 /im maps. However, we conservatively 
establish flux upper limits for COSMOll and for the 
three COSM015 frequencies, which will in turn be im- 
portant to trace the AME flux downturn at lower fre- 
quencies predicted by spinning-dust models. To this aim 
we measure the standard deviation in each map within 
a box of ps 10 square degrees around the Pleiades, ob- 
taining 65, 96, 65 and 72 /iK beam" 1 for 10.9, 12.7, 14.7 
and 16.3 GHz, respectively. After transforming to flux, 
these values are corrected for the flux suppression caused 
by the COSMOSOMAS map-making strategy described 
above. To this aim, we simulate a source with the same 
brightness distribution as the 100 /im DIRBE map in 
the position of the Pleiades, and a flux normalized to 
1 Jy. This source is then degraded to an angular resolu- 
tion of 1.12°, and convolved with the simulated circular 
scanning strategy of COSMOSOMAS, which includes re- 
moval of the first seven harmonics. By measuring the flux 
in this final map, and comparing with the initial flux of 
1 Jy, we conclude that the flux losses are 62% and 66% 
respectively for COSMOll and COSM015. Final upper 
limits are listed in Table [T] at the 99.7% confidence level 
(C.L.). 

3.3.2. WMAP intensity 

The NASA's WMAP satellite have produced full-sky 
maps at frequency bands centered at 22.8, 33.0, 40.7, 60.8 
and 93.5 GHz (K, Ka, Q, V and W), and with angular 
resolutions of 0.85°, 0.65°, 0.51°, 0.35° and 0.2 5°, respec- 
tively. We use the seven-year data products (Jarosik et 



[al~]|201l| , available in the lambda databas^] In Fig- 
ure [3] we represent a composition of the K-band map at 
the position of the Pleiades, with IRAS 100 fim contours 
overplotted, showing an evident spatial correlation be- 
tween both maps. The Pearson correlation coefficient 
calculated between maps at the same angular resolution 
(1°) is r = 0.608. However, the center of the microwave 
emission seems to be shifted ~ 15' to the southwest of 
the position of the stars 23 Tau and 17 Tau, around 
which the bulk of the the infrared emission is seen. Fig- 
ure [5] shows WMAP smoothed maps subtracted from the 
CMB contribution (see next paragraph). AME shows up 
prominently especially in the K band. The bulk of the 



http : //lambda. gsfc .nasa. gov/ 
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Fig. 3. — COSMOSOMAS maps towards the Pleiades reflection nebula, for the four frequencies of this experiment. Note that here we 
are displaying data down to a declination of s« 23.5°. There are data in a different configuration which extends the declination coverage to 
16.7°, but with ~ 20 times less integration time, and therefore with a notably higher noise level. 



emission, which extends ~ 1.5°, seem to be mounted on a 
larger-scale diffuse structure which elongates towards the 
southwest of the maps. Some emission is also hinted at 
in the 0.408-1.4 GHz maps (see Figure [2b. We calculate 
spectral indices between 0.408 and 22.8 ClHz in different 
circular apertures on this diffuse region surrounding the 
Pleiades nebula, obtaining values around —0.5, which are 
compatible with synchrotron emission. However, we can 
not draw a firm conclusion about the possible contribu- 
tion from AME owing to the low flux density and to the 
contamination from the stronger emission arisen within 
the Pleiades dust. Radial profiles of the WMAP K-band 
map at its original angular resolution across the source 
in different directions show the 22.8 GHz emission to ex- 
tend - 1°, being slightly resolved by the 0.85°-FWHM 
beam. This further supports the hypothesis that, even if 
the bulk of the microwave emission may be coming from 
a position ~ 15' southwest of 23 Tau and 17 Tau, there 
may be an important diffuse component originated in a 
wider region encompassing this nebulosity. 

Towards the southwest of the Pleiades stars, at equa- 
torial coordinates (R.A.,Dec.)=55.44° + 21.67°, WMAP 
maps also reveal a source, with a flux of ~ 0.6 Jy at 
22.8 GHz, which has no counterpart in the 1.4 GHz NVSS 
( Condon et al.|1998 ) nor in the WMAP point source cat- 
alogues ( |Gold et al. 2011). It does not have significant 
e mission in th e 0.408 — 1.4 GHz maps nor in the Ha map 
of |Finkbeiner (|2003| ). The flux upper limits that can be 
derived from trie radio data, and the flux value at 22.8 
GHz, are consistent with a free-free spectrum. WMAP 
93.5 GHz data show evidence of thermal dust emission. 
Moreover, as it is shown in Figure [6j this source is visible 
in IRAS 60 and 100 /jm maps as a circular ring with a 
diameter of w 40', as well as in the CO map of |Dame] 
et al. (2001). There are several infrared sources around 
that position listed in the IRAS point source catalogue. 
The spectrum of WMAP and DIRBE data is compatible 
with thermal dust emission. 

For the sake of consistency, we perform the flux extrac- 
tion in maps at the same angular resolution. We then 
use LAMBDA maps degraded to a common resolution of 
1°. Onwards we will consider two separate strategies 
for flux calculation: i) case A: we apply our method for 
flux extraction directly on those maps, which may con- 
tain some CMB contamination; ii) case B: we construct 
CMB-subtracted maps on which our method for flux ex- 
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Fig. 4.— WMAP K-band (22.8 GHz) map with IRAS lOO^m 
data overplotted with contours, in the position of the Pleiades 
nebula. Arrows indicate the positions of the brightest stars in 
the region, and the circle the l°-radius aperture we use for flux 
integration. Temperature units of the WMAP map are indicated 
in the color scale in the top. IRAS contours correspond to intensity 
levels of 20, 30, 50, 100, 200 and 300 MJy sr _1 . A clear correlation 
is seen between the microwave (AME) and the infrared (thermal 
dust emission) intensities. 



traction will be applied. To define the CMB template 
that will be pixel-by-pixel subtracted from the original 
maps, rather than using the full-sky internal linear com- 
bination (ILC) map supplied by the WMAP team, we 
generate our own ILC in a circular area of 30°-radius 
around the Pleiades after masking out point sources by 
applying WMAP's KQ85 mask. We use the same tech- 
nique as the WMAP team, i.e. a weighted linear com- 
bination of the five smoothed frequency maps, in which 
the weights are chosen to minimize the variance of the 
measured temperatures, in order to maintain the CMB 
contribution while minimizing the foreground contribu- 
tion. Pixel temperatures in our ILC are given by Tilc 
= -0.1937 T K -0.0550 T Ka +0.8289 T Q +0.7042 T v 
—0.28437 T\v, and the resulting CMB-subtracted maps 
are shown in Figure [5] 

In order to extract the fluxes we use a dire ct aperture 
integration. This is a well-known technique ( jBennett et 
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Fig. 5. — WMAP and DIRBE maps smoothed to a common resolution of 1°. WMAP maps have been subtracted from the CMB 
contribution given by the 1L(J template (sec text for details). Crosses mark the positions of the same stars shown in Figure |l| The 



l"-radius circular aperture and 1.7-2.0" circular annulus which arc used tor flux extraction arc indicated. Strong AM is evident m'WiVlAb' 
maps, with maximum intensity at 22.8 GHz 
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Fig. 6. — IRAS 100 /im and CO maps of the source 
found towards the southwest of the Pleiades, with coordinates 
(R.A.,Dec.)=55.44° + 21.67°. Crosses mark the positions of the 
infrared sources from the IRAS catalogue. 



al.|1993|IHernandez-Monteagudo fc Rubmo-Martfn|2004| 
Planck Collaborati on et al.|2011[ ), consisting on integrat- 
ing temperatures of all pixels within a given aperture, 
and subtracting a background level which is calculated 
by averaging pixel temperatures enclosed by an external 
annulus. The flux is then given by 



S v = a{v) 



i=l lj - 



En2 rp 



n-2 



(1) 



where n\ and ni are respectively the number of pixels 
in the circular aperture and in the external annulus, and 
Ti and Tj represent the pixel thermodynamic temper- 
atures in the aperture and in the background annulus. 



The function a(u) gives the conversion factor from tem- 
perature to flux, 



a(v) 



2k h Tl 



sinh 



/ hv 



cmb 



(2) 



where h and fcb are the Planck and Boltz mann constants, 

T 



-^cmb 

perature 



Sir 



2.725 K (Mather et al. 
and 



1999) the CMB tem- 



the solid angle subtended by each 
pixel. The associated error bar is calculated through 
the quadratic sum of the variance of each pixel (given by 
WMAP error maps) plus the pixel-to-pixel covariance: 



<r(S v ) = a{v) 



1 



E <w< 



nl 

E« 



mn 2 



E c ^> 



1/2 



(3) 

In this equation 0{ represents the instrumental noise as- 
sociated to each pixel. C(9ij) is the noise correlation 
function, which arises from the smoothing process, and 
is evaluated for any pair of pixels separated an angle 9ij . 
For case A, we also calculate the error associated with the 
CMB primordial fluctuations, by introducing its correla- 
tion function, derived from a theoretical power spectrum 
obtained from the WMAP best-fit cosmological model. 

The radii we use for the aperture and for the inter- 
nal and external circles defining the annulus are shown 
in Figure [5j These are respectively 1.0°, 1.7° and 2.0°, 
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and were chosen from the radial angular profiles of the 
signal. These angular profiles flatten at the distance of 
the background annulus, and the aperture encloses the 
bulk of the microwave and infrared emissions associated 
with the Pleiades. The center coordinates of the aper- 
ture are R.A.= 3 h 44 m 58 s , Dec.= 23°46'48" (J2000). This 
position is chosen to match the maximum emission at 
22.8 GHz, and is located ~ 15' southwest of the mid- 
dle of the line connecting the 17 Tau and 23 Tau stars, 
around which the bulk of the infrared emission is located. 
Final fluxes are shown in Table[T] Error bars include only 
the thermal noise, and the contribution from the CMB 
variance is indicated within brackets for case A. In case 
B the flux error bar includes the ILC noise, which is cal- 
culated in each map pixel through the quadratic sum of 
pixel noises for all bands weighted by the corresponding 
ILC coefficient. The flux at 22.8 GHz calculated in the 
ILC map, which is subtracted in case B, is 0.45±0.11 Jy. 
This value is compatible with 0.51 Jy, which corresponds 
to the standard deviation of the CMB primordial fluctu- 
ations associated with this measurement. 

At 22.8 GHz we get 2.15 ± 0.12 Jy for case B, which is 
two orders of magnitude above the level of free- free emis- 
sion at t his frequency derived from Ha (0.02 Jy; see sec- 
tion 3.1 1. We will see later in section [4jthat this flux is in 
fact dominated by AME, as the contribution from ther- 
mal dust emission is rather low at this frequency. Also, 
considering the upper limit from COSMOll, the spec- 
tral index between 10.9 and 22.8 GHz is ^ii_ 2 3 > 1.12, 
a value consistent with typical spinning dust models. 

3.3.3. WMAP polarization 

Estimates or constraints of the AME polarization can 
help to discriminate between models based on electric 
dipole emission from very small spinning dust grains 
and magnetic dipole radiation from hot ferromagnetic 
grains (see section [l] and references therein) . It is there- 
fore interesting to get estimates of the polarization from 
WMAP maps of Q and U Stokes parameters. Apply- 
ing the sa me aperture photometry technique as in sec- 
tion 3.3.2 in these maps, which show no significant emis- 
sion aFThc position of the Pleiades, we get fluxes of 
Q = 0.076 ± 0.071 Jy and U = -0.056 ± 0.091 Jy 
at 22.8 Ghz. Using these values, together with the 



flux density for case B (see previous section), we de- 
rive a 95% C.L. upper limit on the polarization frac- 
tion, n = y/Q 2 + U 2 /I, of n < 10.8%. This constraint 
is much weaker than n < 1.0% obtained in the much 



brighter region of the Perseus molecular cloud (Lopez- 
Caraballo et al.||2011[ ). 



3.4. Far-infrared data 

Far-infrared data trace the thermal dust emission, 
which can then be extrapolated to microwave fre- 
quencies. We use Zodi-Subtr acted Mission Aver age 



(ZSMA) COBE-DIRBE maps (|Hauser et al 



240 fim (1249 GHz), 140 fim (2141 GHz) and 100 Vm 
(2998 GHz). Although there are DIRBE data at higher 
frequencies, we only consider these three bands, because 
they allow modelling the thermal dust spectrum with 
a single modified black-body curve. Higher-frequency 
bands are dominated by a different grain population and 
their inclusion would make necessary at least two curves. 



We also note that, in spite of IRAS finer angular reso- 
lution, we use here DIRBE data because IRAS lower- 
frequency band is 100 /im, and also because we calculate 
fluxes in a larger angular scale than either of IRAS and 
DIRBE beam sizes. 

In Figure [5] we show the DIRBE maps at the posi- 
tion of the Pleiades smoothed to a resolution of 1°. It is 
again clear from these plots that the microwave emission 
is slightly shifted to the southwest with respect to the in- 
frared emission at 100 fim. This offset is less pronounced 
at 140 fim, and almost disappears at 240 fim. The Pear- 
son correlation coefficients between WMAP 22.8 GHz 
and DIRBE bands from 100 fim to 240 fim are respec- 
tively r = 0.585, 0.722 and 0.821. However, our aper- 
ture encloses the bulk of the flux in both frequency 
ranges. Aperture-photometry fluxes in the three bands 
are quoted in Table [T] We also show the dust-correlated 
emissivities, which are used to quantify the correlation 
between the 100 ^m DIRBE map and WMAP maps. 
These are calculated by a linear fit of the temperature 
(fiK) of WMAP background-subtracted pixels to the 
intensity (MJy sr _1 ) of DIRBE background-subtracted 
pixels using a standard least-squares method. A sig- 
nificant correlation is found, specially between WMAP 
22.8 GHz and DIRBE 100 fim channels, where we get 
4.36 ± 0.17 /xK/(MJy sr _1 ). We note however that this 
dust emissivity is strikingly lower than those obtained 
in dust clouds ( |Davies et al.|[2006| , which usually range 
from ~ 11 to ~ 35 fiK/(MJy sr _1 ), and more similar 
to the 3. 3 ±1.7 ^K/(MJ y sr^ 1 ) found by 



Dickinson et 



al.| (|2007j) in HII regions. |Vidal et al.| ( |201ip got an even 
lower value of 0.2 ±0.1 /xK/(MJy sr _i ) in the translucent 
cloud LDN 1780, where they found evidence of AME. 

4. SPECTRAL ENERGY DISTRIBUTION 

Figures [7] and [8] depict the spectral energy distributions 
(SED) for cases A and B, respectively. In both cases, 
we represent 3-sigma (99.7% C.L.) upper limits at 0.408, 
0.820 and 1.42 GHz, as no signi fican t emission was found 
at these frequencies (see section pT2|) . Upper limits at the 
99.7% C.L. are also shown forCDSMOll and for the 



three channels of COSM015 (see section 3.3.11. 

AME, which neatly shows up in the SEDs as a 
bump at frequencies ^20-60 GHz, is usually well fitted 
by spinning-dust m odels (see section [l] and reference s 
therein). Following Planck Collaboration et al. (2011), 
we consider two complementary spinning-dust compo- 
nents, corresponding respectively to high-density molec- 
ular gas and to low-density atomic gas. Different studies 
of the dynamics of the Pleiades environment agree on 
the existence of these two gas phases, the former be- 
ing probably located within the cluster and the later ly- 
ing ~ 1 pc in front of it (see section [5] and references 
therein). We compute spinning-dust emissivities per hy- 
drogen column density of these two phase s using the SP 
dust. 2 cod<^] ( |Ali-Haimoud et al7 



2009 



Silsbee et al. 



1998 1 at 20111. Rather than fitting all the spinning dust parame 



ters to the observed spectra, we use realistic values and 
fit the amplitude to the spectrum. Assumed physical pa- 
rameters for these two phases, that are used as inputs for 
SPDUST.2, are shown in Table [2j n H is the total hydro- 
gen number density, T g the kinetic gas temperature, x is 

5 http:/ /www. tapir. caltech.edu/~yacine/spdust/spdust. html 
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TABLE 1 

Fluxes and dust-correlated emissivities 







Case A 








Case B 






Flux 


Lvcsiciual tlux 


Oorrcl&tiou 


Flux 


Ivosidual tlux 


OorrclaLion 


(GHz) 


( l\r\ 


{ J y> 


^ivuy sr ; 


(Jy) 


i l\r\ 


/lis./ i^ivuy si ; 


0.408 


< 1.14 


< l.n 


_ 


< 


1.14 


< 1.11 


_ 


0.820 


< 0.89 


< 0.87 


_ 


< 


0.89 


< 0.87 


_ 


1.42 


< 0.51 


< 0.49 




< 


0.51 


< 0.49 


_ 


10.9 


< 1.04 


< 0.87 




< 


0.94 


< 0.91 


_ 


12.7 


< 1.97 


< 1.75 


- 


< 


1.83 


< 1.80 


- 


14.7 


< 1.77 


< 1.48 




< 


1.58 


< 1.56 




16.3 


< 2.43 


< 2.08 




< 


2.20 


< 2.17 




22.8 


2.60±0.06 (± 0.51) 


1.95±0.06 


3.01±0.27 


2.15 


± 0.12 


2.12 ± 0.12 


4.36 ± 0.17 


33.0 


2.55±0.10 (± 1.06) 


1.21±0.12 


0.66±0.17 


1.61 


± 0.15 


1.55 ± 0.15 


2.01 ± 0.09 


40.7 


2.64±0.15 (± 1.59) 


0.64±0.17 


-0.32±0.16 


1.24 


± 0.18 


1.12 ± 0.18 


1.03 ± 0.03 


60.8 


4.71±0.36 (± 3.37) 


0.39±0.40 


-0.77±0.16 


1.75 


± 0.38 


1.23 ± 0.38 


0.59 ± 0.02 


93.5 


9.12±0.89 (± 7.03) 


-0.52±0.97 


-0.25±0.12 


2.94 


± 0.90 


0.37 ± 0.90 


1.10 ± 0.05 


1249.1 


11931±185 


9±394 




11931 ± 185 


-366 ± 195 




2141.4 


23469±249 


-14±595 




23469 ± 249 


618 ± 262 




2997.9 


17959± 89 


1±375 




17959 ± 89 


-47 ± 101 





Notes. Flux upper limits at the 99.7% C.L. level are shown for the 0.408-1.42 GHz and COSMOSOMAS (10.9-16.3 GHz) data. For 
WMAP and DIRBE frequencies fluxes are obtained through an aperture integration, using a radius of 1° for the aperture and subtracting 
a mean background calculated in an annulus between 1.7° and 2.0° around the source. The fluxes on the right (case B) were calculated 
on maps where a pixel-by-pixel subtraction of a CMB internal linear combination map was applied. The error bars include instrumental 
noise, and the numbers within brackets show errors associated with the CMB anisotropics. The residual fluxes represent the level of AME, 
as they were obtained after the subtraction of the rest of the modelled components. We also show the dust-correlated emissivities, which 
indicate the correlation between WMAP and the IRAS 100 jim map. 



the intensity of the radiation field relative to the average 
interstellar radiation field, xu = /% represents the 
hydrogen ionization fraction, xq — nj/nu the ionized 
carbon fractional abundance, y = 2n(H2)/riH the molec- 
ular hydrogen fractional abundance and /3 the average 
dipole moment per atom. Parameters for the grain-size 
distribution are taken fr om line 25 of Table 1 of [Wcin- 
gartner fc Draine| pQOll . 

Different studies of the Pleiades interstellar material 
cited in section[2]provide temperatures aro und T g ~ 20 K 
and d ensities spanning from n ~ 40 cm~ 3 (Ritchey et al. 
2006| to - 400 cm" 3 ( |White||l984"l ). This allows flexi- 
bility when fixing our model parameters. Moreover, we 
note that some of the aforementioned studies focused on 
individual clouds whereas here we study the emission in- 
tegrated over an area encompassing the whole Pleiades 
reflection nebula. The physical parameters we choose for 
the molecular phase a re similar to those o f the mo lecular 
cloud environment of |Draine fc Lazarian (1998b), while 
those of the atomic phase have been adapted to have a 
spectrum peaking at a higher frequency, in order to pro- 
vide a better fit to the 60.8 GHz point. The position of 
the peak of the spinning dust spectrum is very sensitive 
to the intensity of the radiation field. For the molecu- 
lar phase, which is probably associated with the Merope 
molecular cloud, we select a low value of \ = 0.03, in or- 
der to have the peak at ss 23 GHz, matching the WMAP 
K-band flux, which is our maximum measured value. 



Ritchey et al. (12006) have calculated the intensity of 
the radiation field in different lines of sight towards the 
Pleiades; from their values (see their Table 6) in po- 
sitions within our aperture we calculate an average of 
X ~ 16. Using this value and a hydrogen density of 
tin = 30.6 cm~ 3 (also calculated from an average of 
their measurements) , and same values for the rest of the 



SPDUST.2 input parameters, we obtain a spinning dust 
spectrum which peaks at 36.9 GHz and provides a very 
poor fit to the data. If we use for the hydrogen number 
density a value as low as Tin = 0.03 cm~ 3 , while main- 
taining the same radiation field intensity, the spectrum 
peaks at a lower value of 31 GHz, still far from the ob- 
served peak around 22.8 GHz. Therefore, in order to 
reliably fit the data, we have to decrease the intensity of 
the r adiation field f o r the molecular phase. Note however 
that Ritchey et al. ( 2006 ) do not provide measurements 
for the radiation field in any position towards the Merope 
molecular cloud, and therefore it could be the case that 
its intensity is indeed lower. On the other hand, for the 
atomic phase, which probably traces a more diffuse com- 
ponent scattered around the clouds, we adopt x = 10, 
which is closer to the values reported by |Ritchey et al.| 
( 2006 ) , and also provides a better fit to the 60.8 GHz 
point. Once we have obtained the spinning dust spec- 
tra for the molecular and atomic phases from SPDUST.2, 
their amplitudes, which are defined by the hydrogen col- 
umn densities Nn, are fitted to the measured fluxes. 

We model the thermal dust, which is the main contri- 
bution at infrared frequencies, through a modified black- 
body spectrum, I*(T d ) = nor, (j//(2998 GHz))' 3 - 1 B„(T d ), 
with Td being the dust temperature, tioo the optical 
depth at 100 /im and /?d the emissivity spectral in- 
dex. The free-free spectrum is normalized by assuming 
an electron temperature T c = 8000 K, and estimating 
the electron measure, indicated in Table [2j through the 
exti nctio n-corrected Ha emission in the region (see sec- 
tion JO} . This spectrum is represented by the dotted 
lines of Figures [7] and [8] 

In the SED ofcase^ (Figure [7} the CMB is proba- 
bly contaminating our measurements and boosting the 
WMAP V and W-band fluxes. Therefore, in this case 
we also fit to the data a CMB component with an av- 
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Fig. 7. — Top: spectral energy distribution (SED) of the Pleiades 
nebula for case A. Upper limits are shown at the 99.7% C.L. for 
radio and COSMOSOMAS data. WMAP and DIRBE fluxes have 
been obtained through aperture photometry, and are quoted in Ta- 
ble^ The fitted model consists of free-free emission (dotted line), 
thermal dust (dash-dotted line), CMB (dash-triplc-dottcd line) and 
a spinning dust model for a molecular gas phase (dashed line). The 
model parameters are listed in Table [2] Bottom: residual SED af- 
ter subtraction of the free-free and the best-fit CMB and thermal 
dust components. 

erage temperature AT cm b, and neglect the spinning-dust 
corresponding to the atomic gas phase, as it would be 
completely overshadowed by the CMB component. Af- 
ter subtracting the free-free flux corresponding to each 
frequency, we jointly fit, using a non-linear least-squares 
method, our five model parameters, Ah, AT cm b, 7d, tioo 
and /3d, defining the spinning-dust, CMB and thermal 
dust components, to the WMAP and DIRBE fluxes. In 
case B (Figure [8]), we fit the hydrogen column density as- 
sociated with the spinning-dust atomic phase instead of 
AT cm b, as the CMB has been removed directly from the 
map. The fitted parameters and their errors are shown 
in Table [2j and the spectra of all components are rep- 
resented in Figures [7 and [8j We see that these fitted 
spectra are indeed fully compatible with the COSMO- 
SOMAS 3-sigma upper limits. The reduced chi-squares 
are 0.75 and 6.76 for cases A and B, respectively. 

The dust temperature is somewhat higher in case B, 
with the spectrum being steeper in case A. This is a 
consequence of AT cm b in case B being slightly larger 
than the average ILC temperature within the aperture, 
30.2 /iK, which is the effective CMB level that we are sub- 
tracting in case A prior to flux calculation. The fitted 
values for the hydrogen column densities are compatible 
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Fig. 8. — Top: same as in Figure [7] but for case B, i.e. a pixcl- 
by-pixel subtraction of the CMB is performed, and then no CMB 
component is considered in the SED fitting. Fluxes for this case 
are quoted in Table [l] The fitted model consists of free- free emis- 
sion (dotted line), thermal dust (dash-dotted line), and two spin- 
ning dust components: high-density molecular gas (dashed line) 
and low-density atomic gas (dashed-triple-dotted line). The model 
parameters are listed in Table [2] Bottom: residual SED after sub- 
traction of the free-free and best-fit thermal dust components. 

in both cases. These values show that the Pleiades re- 
flection nebula constitute indeed a more diffuse media 
than others were AME have been usually studied, like 
the Perseus and p Ophiuchus molecular clouds, where 



A H 



117 x 10 2U cm" 2 and 171 x 1 2U cm -2 respectively 



(Planck Collaboration et al. 20111. This also becomes 
evident from the reddening mea surem ents in this region, 
Eb-v — 0.1 mag (see section 3.1 1, which are typical 
of diffuse clouds, and lower than the respective values 
in the Perseus and p Ophiuchus molecular clouds, and 
even lower than the reddening of the translucent c loud 
LDN1780, £ B -v - 0.6, where iVidal et all poTTb re- 



cently found AME. Our fitted values for the hydrogen 
column densities are also of the same order of those de- 



rived by applying the scaling relation of |Bohlin et al.| 



fl978j , (A H + N K2 )/E B -v = 5.8 x 10 2i cm^ mag" 1 , 
and taking for the reddening E-&— v = 0.1 mag. It may 
also be noted that by dividing the fitted hydrogen column 
density of the molecular phase by the hydrogen number 
density, we get an estimate for the depth of the cloud 
of ~ 0.7 pc, a value that agrees with the characteristic 
length scale of this cloud. 
It is also interesting to co mpare our results wi th those 



presented in Figure 13 of Vidal et al. (2011), where 



they plot emissivity at 31 GHz (intensity at 31 GHz 
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TABLE 2 




Model parameters 


Te (K) 


8000 


EM [era 


~ 6 pc) 0.267 


Molecular Atomic 



nn (cm 3 ) 


300 200 


T s (K) 


20 1000 


X 


0.03 10 


ih (ppm) 


9.2 373 


x c (ppm) 


1 100 


y 


1 0.1 


/3 (D) 


9.34 9.34 


Case A 


JV H (10 20 cm" 2 ) 


6.94 ± 0.22 


T100 


(6.09 ± 0.06)xl0~ 4 


/3d 


2.29 ± 0.02 


Td (K) 


20.12 ± 0.03 


AT cmb (fj,K) 


42.2 ± 1.9 


Case B 


N H (10 20 cm -2 ) 


6.60 ± 0.11 0.30 ± 0.01 


TlOO 


(3.302 ± 0.004) xl0~ 4 


/3d 


1.869 ± 0.004 


^d (K) 


22.008 ± 0.005 



Notes. The two values on the top of the table define the spectrum 
of the free-free emission. T c is an assumed value, and EM is cal- 
culated fro m th e Ho intensity corrected from Galactic extinction 
(see section |3.1| . Parameters from njj to j3 (see section [4] for an 
explanation ol their meaning) are used as inputs for the spinning- 
dust models. They are fixed at typical values for the two different 
gas phases (molecular and atomic). The rest of the parameters 
are jointly fitted to the spectrum. Nn defines the amplitude of 
the spinning dust component; T100, /3d and define the thermal 
dust spectrum; and AT cm t> ls the average value of the CMB tem- 
perature within the aperture. This CMB component is fitted only 
in case A, as in case B we calculate the fluxes on CMB-subtracted 
maps. Also, in case A no atomic gas phase producing spinning-dust 
emission was considered. 

divided by hydrogen column density) versus hydrogen 
column density. Dividing the total flux at 33 GHz 
for case B by the solid angle of the aperture and by 
the corresponding fitted column density we get £33 = 
(2.55 ± 0.28) x 10~ 24 MJy sr" 1 cm 2 . In order to scale 
to 31 GHz we use our fitted model, and obtain £31 = 
(3.03 ±0.33) x 10~ 24 MJy sr" 1 cm 2 . This value roughly 
follows the trend shown in that figure, but falls below 
the fitted line, which for ATjj = 6.60 x 10 20 cm~ 2 predicts 
~ 5.7 x 10 -24 MJy sr _1 cm 2 . It therefore seems that the 
Pleiades nebula, in addition to having a notable lower 
hydrogen column density to the regions shown in that 
figure, and actually to the majority of the regions where 
AME has been previously detected, also has a lower emis- 
sivity per hydrogen nucleon. 

The residual fluxes, obtained after subtracting all mod- 
elled components except the spinning-dust are quoted in 
Table [I] and plotted in bottom panels of Figures [7] and 
[HJ According to our fitted model for case B, the ther- 
mal dust contributes at 22.8 GHz with a flux of 0.012 Jy, 
whereas the level of the free-free at this frequency was 
0.019 Jy. Then, the residual flux at this frequency, which 
can be attributed to AME, is 2.12±0.12 Jy (17.7a; mod- 
elling errors included). 

5. CONCLUSIONS 



We have presented evidence for anomalous microwave 
emission in the position of the nearby (125 pc) Pleiades 
star cluster, one of the most studied regions in the sky. 
This AME is originated in the dust grains that make up 
the Pleiades reflection nebula, 



whose thermal infrared 
Us- 



emission was discovered by Castelaz et al. (19871 
ing data from the seven-year release of WIVlAP, here we 
have measured a total flux density towards this nebulos- 
ity of 2.15 ± 0.12 Jy at 22.8 GHz by applying an aper- 
ture integration in a l°-radius circle around the position 
(R.A.,Dec.)=56.24° + 23.78° (J2000). COSMOSOMAS 
maps at 11-17 GHz do not show evidence for emission 
at this position. Nevertheless, we have set upper lim- 
its (99.7% C.L.) of 0.94 Jy at 10.9 GHz and 1.58 Jy at 
13.6 GHz, which help to trace the downturn of the AME 
spectrum at these frequencies predicted by models based 
on electric dipole emission from fast-spinning dust grains. 

In order to quantify the total AME flux we have con- 
sidered other emission mechanisms. We have estimated 
a free-free flux of 0.019 Jy, at the same re ference fre 
que ncy of 22.8 GHz, from the Ha map of Schlegel et 



al. ( 1998 ) , which has been previously corrected from 
Galactic extinction by using reddening values measured 
towards different stars. Infrared fluxes obtained from 
the COBE-DIRBE experiment have been used to trace 
the thermal dust emission, from which we extrapolate 
to microwave frequencies and obtain a flux of 0.012 Jy. 
Therefore, both the free-free and thermal dust emissions 
have little impact in our measurement. After deducting 
these two components we have obtained a residual flux of 
2.12 ± 0.12 Jy, which is a detection of AME at the 17.7<r 
level. 

The signal at 22.8 GHz is slightly resolved by WMAP 
K-band beam (0.85° FWHM). This indicates that, al- 
though the bulk of the emission is probably coming from 
the region near the stars 23 Tau and 17 Tau, where 
the strongest infrared emission is also detected, there 
is probably an important diffuse component integrated 
over the whole nebulosity. Observations with microwave 
experiments at higher angular resolution are crucial to 
elucidate the relative contribution to the total observed 
flux of the different individual clouds that constitute the 
Pleiades reflection nebula. They may also help to better 
trace the correlation with the thermal dust, and to shed 
light on the properties of the grain population that may 
be originating this AME. 

After deducting the extrapolated free-free fluxes, we 
have fitted a joint model consisting of a modified black- 
body spectrum, tracing the thermal dust emission, and 
a spinning-dust spectrum which traces the AME, to the 
WMAP and DIRBE fluxes. For the spinning-dust com- 
ponent we have considered two complementary phases, 
one dominated by high-density molecular gas and the 
other by lower-density atomic gas. Their spinning-dust 
spectra have been obtained with SPDUST.2, by using re- 
alistic values for the physical parameters that are inputs 
to this code. We then fitted the amplitude of the spec- 
tra to the data, together with the three parameters that 
define the thermal dust model. The resulting spinning- 
dust fit successfully traces the WMAP and DIRBE data, 
while being compatible with the COSMOSOMAS upper 
limits. 

Microwave temperatures at 22.8 GHz are clearly cor- 
related with 100 /jm infrared intensities, the dust emis- 
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sivity being 4.36 ± 0.17 /zK/(MJy sr^ 1 ). This value is 
lower than what is typical in cool dust clouds, and more 
characteristic of HII regions. The physical properties of 
the Pleiades nebula show that this is a much less opaque 
medium than others where AME has been generally stud- 
ied, the reason for this detection likely being its proxim- 
ity. For instance, the optical depth and the hydrogen 
column density of the Perseus mole cular cloud are more 



than an order of magnitude highe r (Watson et al. 2005 
Planck Collaboration et al. 2011), which explains why 



the AME flux at 22.8 GHz is there twenty times larger 
than in the Pleiades despite being at twice their distance. 
This potentially makes the Pleiades an attractive envi- 
ronment to study the AME emission in different physical 
conditions than those generally explored up to now. 
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where we got the 0.820 GHz data. We acknowledge the 
use of the Legacy Archive for Microwave Background 
Data Analysis (LAMBDA). Support for LAMBDA is 
provided by the NASA Office of Space Science. Some of 
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